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Abstract 

We analyze the high temperature (or classical) limit of the Casimir effect. A useful 
quantity which arises naturally in our discussion is the "relative Casimir energy", which we de- 
fine for a configuration of disjoint conducting boundaries of arbitrary shapes, as the difference of 
Casimir energies between the given configuration and a configuration with the same boundaries 
infinitely far apart. Using path integration techniques, we show that the relative Casimir energy 
vanishes exponentially fast in temperature. This is consistent with a simple physical argument 
based on Kirchhoff's law. As a result the "relative Casimir entropy", which we define in an 
obviously analogous manner, tends, in the classical limit, to a finite asymptotic value which 
depends only on the geometry of the boundaries. Thus the Casimir force between disjoint pieces 
of the boundary, in the classical limit, is entropy driven and is governed by a dimensionless 
number characterizing the geometry of the cavity. Contributions to the Casimir thermodynam- 
ical quantities due to each individual connected component of the boundary exhibit logarithmic 
deviations in temperature from the behavior just described. These logarithmic deviations seem 
to arise due to our difficulty to separate the Casimir energy (and the other thermodynamical 
quantities) from the "electromagnetic" self-energy of each of the connected components of the 
boundary in a well defined manner. Our approach to the Casimir effect is not to impose sharp 
boundary conditions on the fluctuating field, but rather take into consideration its interaction 
with the plasma of "charge carriers" in the boundary, with the plasma frequency playing the 
role of a physical UV cutoff. This also allows us to analyze deviations from a perfect conductor 
behavior. 
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1 Introduction 



There has been considerable interest in the Casimir effect [l] lately 0]. Recent discus- 
sions of this effect and its variants range over several areas of physics [[| [|, [7], [^, 
|H], PH| . For recent books and reviews see |L2|, [13|, [14], [15], ^6| . For earlier accounts 
see 0, 0|. 

The experimental status of the effect was tenuous until 1997 when vindication of 
the theoretical prediction to an accuracy of a few percent was reported 

Over the years, various methods have been applied to analyze the Casimir effect. 



To name a few, we mention the "mode summation technique" |L2], [L3|, [L7|, [18] (which 
was the method originally used by Casimir), analyses of Green's functions and of the 
energy-momentum tensor in the presence of boundaries | 



25, 



and 



zeta function techniques |27|, Effil . These techniques possess flexibility which permits, 
among other things, the handling of temperature dependence (in vacuum and in 
dielectrics) 12T 



various 



% m M, BO, 31, 32, gl, of fields with finite mass [31, 35 



geometrical setups || ^ [7J, [22], |24|, |25], |33|, |36], , and a better insight into the 



Casimir force (macroscopic) in relation to microscopic "long range forces" [|14 



Finally, we mention here the path integral approach H |6|, |7|, ||, [38[ which seems 
to be well suited to handle systems with boundaries of arbitrary shape as well as 
systems with moving boundaries . 

An important aspect of the Casimir effect is its behavior at finite temperatures. 



Finite temperature studies were carried out by several authors, e.g., p(| [21|, |26|, ^j, 
50 , |3T] , p4]| . In particular, we mention in this context the thorough analysis by Balian 
and Duplantier |B2| , who used multiple scattering expansions to study the Casimir 
effect at finite temperature in the presence of conducting boundaries of arbitrary 



smooth shapes. One of the several issues studied in [3!| was the classical (i.e., high 



temperature) limit of the effect (see Section 6 in f32 



In this paper we address the classical (or high temperature) limit of the Casimir 
effect due to boundaries of arbitrary shapes using path integration. Our results are 
in agreement with those of E|. The classical limit is attained upon subtraction of 
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the "self-energies" of the boundaries. The need of this subtraction for obtaining the 
classical limit of the Casimir effect leasds us to define in (Section 2 of ) this paper 
what we call the "relative Casimir free-energy", the "relative Casimir energy" and 
the "relative Casimir entropy". In fact, based on Kirchhoff's law, we can explain the 
classical limit of the Casimir effect intuitively in terms of these newly defined relative 
quantities. 



To get oriented, let us briefly review (hopefully in a self contained manner) the 
standard Casimir problem, viz., radiation field confined between two parallel conduct- 
ing plates, at zero and at finite temperature. We will write down explicit expressions 
for the Casimir energy, free energy and entropy. These are displayed in Fig. 1, which 
gives the behavior of these quantities for all temperatures. At the end of this sec- 
tion we will discuss the classical limit of the Casimir effect in this simple geometry. 
(Readers familiar with these standard aspects of the Casimir effect may skip directly 
to the next section, but are urged to take a look in Fig. 1 on the way.) In Appendix 
C we present an alternative treatment of this problem, based on the path integral 
formalism, following closely @, 0- HI • 

For simplicity, we will employ, for the rest of this section, the so called "mode 
summation technique". (For more details see, e.g., [[[(], |12|, [13|, [17], [TJ|. Here we follow 



closely flllj.) We consider the radiation field (in 3 + 1 space-time dimensions) confined 

between two conducting parallel plates. (In later sections, where we discuss the high 

temperature limit of the Casimir effect for cavities of arbitrary shapes, we will limit 

ourselves to massless scalar fields[] for simplicity.) The sizes of the plates' edges are 

L. The parallel plates coincide with the planes at z = and z = d. We are interested 

of course in the case where L » d, and in the end will take the limit L — > oo. 

The energy tied down in zero point fluctuations in the mode k is, in obvious 

1 This simplification does not spoil the agreement between our results and the results in Section 
6 of p2] which discusses the electromagnetic case. 
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notation (k includes the polarization and the zero in the argument relates to T = 0), 

E k (0) = hu k = hc\k\ . (1) 

The total energy density of the EM field with the conducting plates as boundaries 
is (k± is the magnitude of the wave vector perpendicular to the z axis, i.e., parallel 
to the plates p 



£(d,T = 0) 



he r°° 



L 2 d 2nd Jo 



k\dk\ 



+ 5Z \k r , 



A i 



1,2 | 1 2 i 2 i 2 , m2?I " 2 

K x + K y , K m — K ± H : 



m=l 



m = 1, 2, ... 



(2) 



The energy density, in dimensionless units, e(0), is given by |TO 



£{d,0) 



7T 



2 ^C 



2d 4 



e(0) = De(0); 



e(0) 



X 



+ v^T 



m=l 



(3) 



(4) 



Equation (H) defines the energy density scale D. The energy density due to vacuum 
fluctuations of the radiation field in an arbitrarily large volume V = L 3 - which serves 
as the reference (unconstrained) system - is 



£(oo, 0) = — 2 \ huJ k = D£oo{0), 

k 

oo oo 

£oo(0) = J xdx J dm\/ x 1 + m 2 . 
o o 

The Casimir energy density, at T = 0, in dimensionless units, is given by 

£ c (0) =e(0)- eoo (0). 



(5) 



(6) 



Both e(0) and £oo(0) diverge. e c is commonly |12|, evaluated by a physically 
justifiable regularization technique. Thus a wave vector dependent function, r(k/k c ), 
is introduced into the above integrals such that for k » k c r —>■ while r — > 1 for 
k << k c , thereby rendering the integrals convergent. A simple choice is (a = l/k c ), 



r(ak) = exp [— ak] 



(7) 



a is allowed to go to zero at the end of the calculations - for all the terms together - 
thus the sum is "regularized". The details of the calculations will not be given here 
(cf. PH| , |T7| [T8|j.) In Appendix C, where we present an alternative analysis of this 
problem, we employ a different regularization scheme for the Casimir energy, which 
seems to us more natural physically. 
The result for our case is 



e e {0) 



(27T) 



rC(4) 



or, in physical units, 



S c (d, 0) = De c (Q) 



360 



7T 2 TlC 

~720d 4 



(9) 



The force f c (d, 0) per unit area between the plates is now calculable from the regu- 
larized Casimir energy (H), yielding the well-known Casimir force |]] 

1 d{L 2 dS c ) 7t 2 hc 



fc(d,0) 



(10) 



L 2 dd 240(i 4 ' 

The finite temperature problem is quite similar. Now the zero point energy is 
supplemented by the thermal energy. Thus the (average) energy tied down in the 
k-th mode is 



E k {T) = {l/2)hu k 



huok 



exp(/3hujk) — 1 
Correspondingly, the energy density is 

he 



: coth 



£(d,T) 



2nd 



k i dk i 



\k±\ f /3tuj(k A 
— — coth ( — 



+ |^m|coth( 



j3huj(k r , 



m=l 



(12) 



Returning to our dimensionless units we may write the energy density as a sum of a 
zero temperature part plus a temperature dependent part: 



e{d,T) =e(0)+u(d,T), 
/(0) 



u(d, T) 



+ E/( 

m=l 



m) 



(13) 
(14) 



f(m) = dy y 2 n(y,T), 



with 

k B T c = Tic- . (16) 

Observe that for dimensional reasons, e(d,T),u(d,T) and n(y,T) (as well as all the 
other thermodynamical quantities which we will define below) depend on d and T 
through the combination Td oc T/T c only. Thus, in particular, e(d,T = 0) = e(0) is 
independent of d. 

We see from ( |T5D that T c plays the role of a characteristic (or lower cutoff) tem- 
perature of the geometrical setting, separating high and low temperatures. The cor- 
responding expression for the unconstrained system is 

£■(00, T) = e(0) + u(oo,T), (17) 

with 

oo 

m(oo,T) = / dm f(m) . 



Evaluating the sum in flH| ) via the Poisson summation formula [HJ gives (/i = 

27T777.) 

OO 

e c (t) = e(d, T) - e(oo, T) = -At 3 ^ -coth(fy/)csch 2 (£//), t = vr— . (18) 

m=l A* "'c 

As a consistency check we observe that for t << 1, e c — > s c (0), the expression 
in (§). The finite temperature Casimir energy £ c (i), i.e., the difference between the 
expectation values of the Hamiltonians of the constrained and unconstrained systems, 
is displayed in Fig.l. 

The Casimir free energy is calculated as follows: The partition function (for one 
mode) is 

Z k = ± exp [-(3hu k {n + 1/2)] = ^tS^L (19) 
t^x 1 - ex P {-priu k \ 

and thus the free energy for the k-th mode is 

F k = ~k B T log Z k = ^hu k + k B T log(l - exp [-(3TiL0 k \) . (20) 



Hence the expression for the free energy density of the constrained system is 

he 



F(d,T)=£(d,0) + k B T^J k A 



dk 



log(l - e -^ cfc± ) + 2J2 Ml - e' phckm ] 

m=l 

(21) 

A corresponding equation holds for the unconstrained system, thereby leading to the 
dimensionless expression for the Casimir free energy density, C , given by 

F c {T) = D<f> c (t), (22) 

where D was defined in (|3|), and 

(f) c (t) = -2t -3 fcoth(t/x) + (t/i)csch 2 (t/i)j . (23) 

m=l A* 

The Casimir free energy, (j) c {t), is displayed in Fig. 1. As a check, note that for 
t « 1, 4> c {t) — > e c (0), as it should. 



At this juncture it is natural to consider Casimir's entropy |10| , [2T1 . This entropy, 
cr c (t), is defined in our dimensionless units by 

<t>c{t)=e c {t)-tcj c {t). (24) 

—<j c {t) is displayed in Fig. 1. 

At high temperatures, t » 1, e c (£) falls off exponentially with t, a c (t) approaches 
a constant value (independent of t and, of course, of h) while <f> c (t) becomes propor- 
tional to t. (An explicit calculation of this behavior is deferred to Appendix A.) Re- 
expressing these results in standard physical dimensions we get, at this limit (t » 1) 

S c ^0, F C ^-TS C , S c ^^&. (25) 

2'% GT 

Thus the entropy is proportional to the area of the plates scaled by d 2 . In view of 
this fact we conclude that in the classical limit (high temperatures) the Casimir force 
is purely entropic. 
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Figure 1: Casimir's Energy e c (a), Free Energy C (b), and Entropy —o c (c) as a 
function of temperature, t. Note that graphs (a) and (b) coincide at energy —0.0028 ~ 
— 1/360, consistent with (H). 

It is clear from Fig. 1 that the Casimir energy vanishes at high temperatures (and 
according to (|18|) it vanishes exponentially fast). This, we contend, characterizes the 
classical limit where the energy equipartition "theorem" and the validity of the one 
to one correspondence between the states of the constrained system and the states 
of the unconstrained system imply this result. This is discussed in more detail in 
Section 2. 

This concludes the general introductory review. The rest of the paper is organized 
as follows. Section 2 includes the central point of the paper, viz., an exposition 
of the classical limit and a study of its implications. It is based on the argument 
that the known [0] one to one correspondence between states of constrained systems 
and essentially "free" ones implies the vanishing of the Casimir energy. This is a 
somewhat subtle issue. To be precise, what vanishes in the classical limit is what we 
define as the relative Casimir energy, namely, the Casimir energy associated with a 
given disconnected set of boundaries, from which we subtract the Casimir energy for 
a system in which all these boundaries are infinitely separated from one another. In 
particular, this statement does not hold for a boundary made of a single connected 
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piece. In contrast, the relative Casimir entropy does not vanish in the classical limit. 
We also summarize in Section 2 the high temperature results of Section 6 of |32| . 
which support these assertions about the classical limit for systems constrained by 
conducting boundaries of arbitrary smooth shapes. In Section 3 we study, in the 
path integral formalism, the Casimir effect associated with a massless scalar field 
(f>(x) (as a simple analog of the electromagnetic field), confined in a static reflecting 
cavity of arbitrary shape. Section 3 is divided into several subsections. In Section 
3.1 we study the zero temperature effect. First, we review the formalism of |5], |7], |J 
assuming sharp boundary conditions on all frequency modes (Section 3.1.1). These 
boundary conditions on the field modes are of course an approximate macroscopic 
manifestation of the interaction of the "plasma" of "charge carriers" that "live" in 
the conducting boundary, with the field <ft, in the limit that the boundary is a perfect 
conductor. Thus, in Section 3.1.2 we adopt a more realistic approach and take into 
account the "susceptibility" of the boundary medium, with its plasma frequency uj p 
playing the role of a physical UV cutoff. In the limit of infinite u p we recover, of 
course, the effect of sharp boundary conditions. However, by keeping uj p large but 
finite, we can also calculate corrections to the Casimir effect due to deviations from 
perfect conductor behavior. In Section 3.2 we discuss the finite temperature case. 
The relevant temperature range is ksT « hu p . In Section 3.2.1 we study the 
classical limit /cbT c < k^T << hu p , where T c is a lower (IR) cutoff temperature 
analogous to fll6l). For this model, we are able to verify the arguments of Section 2, 
that the classical limit for this system results in a vanishing (relative) Casimir energy 
and a geometry dependent (relative) Casimir entropy. This we can do only when the 
boundaries are made of disjoint pieces. We find that contributions to the Casimir free 
energy and entropy due to each connected component of the boundary individually 
exhibit logarithmic deviations in temperature from the behavior just described. In 
the conclusions which are given in Section 4 we review the argument (cf. |12|) for the 
existence of zero point fluctuations as viewed from our vantage point. 

Some technical details are relegated to the appendices. In Appendix A we provide 
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an explicit calculation of the exponential fall-off of e c (t) seen in Fig.l (which we 
mentioned following Eq. ( p4[) ). We show that a naive high temperature expansion 
for the Casimir energy in powers of 1/T does not have any corrections to the classical 
value. In Appendix B we verify the consistency of uj p (T = 0) as a UV cutoff at 
high temperatures. Appendix C provides a derivation of the Casimir effect for the 
canonical parallel plate geometry via the path integral formalism of Section 3. We 
obtain all the results of Section 1 at zero and finite temperature, as well as the classical 
limit. Appendix D contains some technical details concerning bounds on terms in a 
certain expansion we carry out in Section 3. 
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2 The Classical Limit - Discussion 



An intuitive understanding of the vanishing of the Casimir energy of the parallel plate 
geometry for T > T c ( cf. Fig. 1) may be gained via the following reasoning. Loosely 
speaking || the total number of normal modes (per unit volume) in our confined 
system is unchanged upon changing d (in the case under study in Section Thus 
moving the walls adiabatically leads to shifts in the levels - not to appearance (or 
disappearance) of levels (modes), and there is an (intuitive) one to one correspondence 
between the levels regardless of the size of d. This, coupled with the classically valid 
equipartitioning of energy, which of course means that each harmonic oscillator-like 
mode holds k^T amount of energy, implies that the (bulk) energy density (in this 
classical limit) is unchanged upon changing d. 

Thus, in the classical limit, defined by the validity of the Rayleigh Jeans (RJ) law 
or the equipartioning of the energy j|0| , the difference in the energy density between 
the constrained (d « L) and the unconstrained (d ~ L) cases, namely, the Casimir 
energy, is nil. 

Thus we have, in conformity with the result of the exact calculations as exhibited 
in Fig. 1, that, 

lim e c (d,T) = 0. (26) 

J-c 

These considerations are readily generalizable to systems made of several disjoint 
conducting boundaries of arbitrary shapes (the parallel plate geometry being an ex- 
ample of such a system). In the next section we will prove, by an explicit calculation, 
that the relative Casimir energy £™ 1 , which we define as the difference (in obvious 
notation) 

S r c el (T) = E c (T)-E™(T) (27) 

between the Casimir energy associated with a given disconnected set of boundaries and 

the Casimir energy for a system in which all these boundaries are infinitely separated 

2 The calculations of Appendix A constitute an explicit 'proof of this claim (within the regular- 
ization scheme employed there). 
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from one another (but are otherwise unchanged), vanishes in the classical limit 

hm^(T) = 0, (2* 



T 

t£ ►OO 

J-c 



where T c is a characteristic temperature of the constrained system, analogous to (|X6|) . 
(We tacitly assume that T c is the maximal parameter of this kind, if several such 
parameters are needed to describe the constrained system.) 

In discussions of black body radiation, Kirchhoff 's law (that the ratio of emissiv- 
ity to absorptivity of all bodies is a universal function of the wavelengh A and the 
temperature T) is used [0 to infer that U, the total electromagnetic energy density 
in a cavity at thermal equilibrium, is a function of T only, 

U = U(T) . (29) 

In particular, it means that U (T) is independent of the shape of the cavity. It is clear 
that (|29|) is true only for T > T c . (E.g., at T = it is, strictly speaking, never true.) 
For this reason we refer to ( p8|) as "Kirchhoff 's theorem" 

Coming back to the parallel plate geometry, it follows from (28), via (^5|), that 



the Casimir force f c = —dT c jdd between the plates is purely entropic in the classical 
limit. Furthermore, in this limit, the Casimir entropy S c (defined as the difference in 
entropies between the constrained and free systems) is independent of temperature 



and is, hence, purely geometric [10 



More generally, as a result of "Kirchhoff 's theorem" (|2~8"D , the Casimir forces be- 
tween disconnected boundaries become purely entropic in the classical limit. Con- 
sequently, the relative Casimir entopy S^ el (T), defined (in analogy with the parallel 
plate case) as the difference 

S r c e \T) = S C (T) - S C °°(T) , (30) 

between the Casimir entropy S c associated with a given disconnected set of bound- 
aries, and the Casimir entropy of a system in which all these boundaries are 

3 That U(T) depends on the shape of the cavity at low temperature, which thus sets limitations 
on the validity of Kirchhoff 's law, is of course well known. See, e.g., the comment made on p. 78 in 
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infinitely separated from one another (but are otherwise unchanged), becomes, in the 
classical limit, independent of temperature and hence, purely geometric. 

This is consistent, of course, with physical intuition, because at temperatures 
T » T c the dominant contribution to entropy is from short wavelength modes 
(A ~ 1/T). These are insensitive to the boundaries when the dimension of the cavity 



much exceeds these wavelengths |42|. This dominant contribution to the total en- 
tropy cancels out when one calculates the Casimir entropy, namely, the difference in 
entropies between the unconstrained and constrained systems. The leading contribu- 
tion to the Casimir entropy, then, is essentially determined by the long wavelengths 
which are geometry dependent, i.e., which relate to the "shape" of the cavity and are 
independent of temperature. 

Let us formalize this argument somewhat further: The entropy of a non-interacting 
Bose gas, such as the gas of photons confined in the cavity, is 



S(T) 



1 

T 



dui p{ui) 



- k B T log ( 1 - e 



(31) 



where p(uj) is the (extensive) density of modes. Thus, for photons (inside a very large 
volume V) pf ree (u) = Vu 2 /n 2 c 3 , whereas for radiation modes between the parallel 
plates of the previous section 



C7T 



/ud)p fT 



1/2 + e((iod/cn) - 



n 



n=l 



and so on. 

Based on ( |5T| ) we can express the Casimir entropy as 

Two 



S C {T) = - J duj (p(u) - P/ P6e (a;)) 



k B T log (l 



-/3hw 



(32) 



Clearly, the density of modes p(uj) in the cavity differs significantly from the den- 
sity of modes pf ree {u) of the free system only at long wavelengths, i.e., they differ 
significantly only for frequencies below some IR cutoff frequency uj c which is deter- 
mined by the shape of the cavity. Thus, the spectral integral in (|32|) is dominated by 
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it) < uj c . This allows us to define the high temperature (classical) limit as such that 
(5hui < (3Tiuj c = T c /T « 1. In this limit (and due to the fact that p — Pf ree cuts the 
integral in (|3"2]) above uj c ) we see that 

oo 

lim S c = k B \ du (p(uj) - p free (uj)) [1 - log (f3hu)] . (33) 

We will treat 

5p(u) = p(uj) - Pfree(uj) (34) 

as an effective Casimir density of states, with normalization^ 

oo 

C = J ' duj 5p(uj) . (35) 



Following Section 6 of |32[1 , we define an average frequency u through 

1 °° 

log u = — J du 5p(uj) log uj . (36) 


In terms of these quantities, we can write (^) as 

lim S c = C (1 - log (J3hQ)) . (37) 

I c 

Balian and Duplantier were able to show that the normalization C associated with 
electromagnetic fluctuations is 

C=— / d 2 a(^ + + -n, (38) 

128 7T J boundaries \ Ki K2 ix\tX2 ) 

where R\ and R2 are the local curvature radii of the boundary, and n is the genus of 
the boundary. Using the Gauss-Bonnet formula 

1 f d 2 a , , 

n = n c / , (39) 

Air J boundaries R\R,2 



4 As we already discussed, there is a one to one correspondence between the modes of the con- 
strained and the free systems. However, we cannot use that correspondence to break (|3^), which is 
supposedly a convergent integral, into the difference of separate integrals over p and p/ ree , which 
are equal (due to the correspondence between states), but infinite. 
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where n c is the number of connected pieces of the boundary, we can write ( |38| ) alter- 
natively as 

1 f ,2 ( 3 3 34 \ , An . 

C = — — / d 2 a —J + -g + — — - n c . (40) 



1287T . /boundaries V-^l R\ R1R2 

An obvious feature of fl3"8| ) and (f40|) is the absence of a term proportional to the 
area 

A= f d 2 a. 

J boundaries 

This is peculiar to electromagnetic fluctuations, due to cancellation between electric 
and magnetic modes.f] In the case of fluctuations of a relativistic massless scalar field, 
for example, which we discuss in Section 3, Balian and Blochf] showed that 

WW ) = p(.) - P /r «M = =F £ + ^ «Pa + ^) + O (i) , (41) 

where the '— ' sign occurs for Dirichlet boundary conditions and the '+' sign, for Neu- 
mann boundary conditions.[] Thus, the normalization C for scalar field fluctuations, 
will have a (divergent, due to the uj integration) area contribution. 

From these considerations we see that the Casimir entropy in the classical limit, 
as given by (0), depends on the geometry of the boundaries, but it also depends 
logarithmically on temperature. It is not just a geometric dependent constant as 
required by "Kirchhoff's theorm". The Casimir entropy (^) can be derived in the 
usual manner from the Casimir free energy 

^--CTlog^V (42) 



which is precisely the high temperature expression obtained in [32] via the mutiple 
scattering expansion. We can interpret fl4"2"|) as the high temperature limit of the free 
energy of C non-interacting harmonic oscillators of frequency lD. From ( |42l) we thus 
find that the Casimir energy is 

£c= d jm = CT, (43) 



5 For more details, see Section 4 (and in particular, Sections 4.2 and 4.3) in p2y . 

6 See Eq.(1.7) in the first reference cited in [F33|. Our p(uj) is 2w times the expression in j33|. 

7 The fact that Dirichlet and Neumann boundary conditions induce opposite signs for the area 
term in (^l]) is essentially the reason why, in the case of electromagnetic fluctuations, the area term 
contributions of electric and magnetic modes cancel each other. 
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and it does not vanish. It is of course a statement of classical equipartition: according 
to fl35|) we can interpret C as the number of additional modes of finite frequency 
created by introducing the conducting boundaries, each of which carries thermal 
energy k B T. 

How do we reconcile these high temperature behaviors of S c , S c and T c with our 
expectations of the classical limit? It is at this point that the relative thermodynamic 
quantities which we defined above become useful: It is clear from the explicit formulas 



(|38| ) and ([H]) that C depends only on the geometrical features of each individual con- 
nected piece of the boundary, as well as on the overall topology of their arrangement 
in space, but not on the relative distances and orientations of these components. (On 



the other hand, it is clear from j32| that the average frequency to, does depend on the 
the relative distances and orientations of these components; otherwise, there will be 
no Casimir forces between them.) Thus, we immediately conclude that the relative 



Casimir entropy (|30|) tends in the classical limit to 



Jim S: el = C\og(^\ (44) 



J-c x 



(where lj(oo) is the average frequency (^) of the system with all boundaries far 
apart), which indeed depends only on the geometry of the boundaries. Similarly, 
the relative Casimir energy (|3), S r c el (T) = S C (T) - £ C °°(T), vanishes in the classical 
limit, because the constrained system in which all the same disjoint boundaries are 
infinitely far apart, has the same coefficient C. "Kirchhoff's theorem" (|28| ) holds. 
It follows then that the relative Casimir free energy T r c el \ namely, the difference 

=T C - (45) 

between the Casimir free energy of the constrained system and a system with the 
same boundaries infinitely far apart, is proportional to the temperature in the classical 
limit, 

T rel = _ TS rel _ ( 4 g) 

Hence, the Casimir forces between disjoint boundaries are purely entropic in the 
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classical limit, as was mentioned above. This reasoning is born out by the explicit 
calculations presented in the next section and in Appendices A and C. 

An interesting point to note is that sign of C can be either positive or negative^, as 
should be clear from the explicit expressions ([38]) and (f40)|32l. For example, a single 



toral boundary will have a positive C if the hole in the center of the torus is very 
small (i.e., if the torus is a very fat tube), or if the torus is very narrow compared 
to the radius of the hole (i.e., the torus is a very thin tube). A torus inbetween 
these domains, will have a negative C. (For more details see [0.) This means 
that the classical relative entropy S^. el ([44]) (being a difference of entropies) need not 
be positive, and, also, that the Casimir forces between the disjoint boundaries may 
be either attractive or repulsive, depending on the geometry (and topology) of the 
boundary pieces. (The signs of S r c el and of the forces depend of course also on u.) 

We are now in a position to interpret the zero point energy (zpe) as a contribution 
necessary to assure that at high temperatures , the energy, U, is a function of T only 
as, indeed, was noted long ago (in 1913) by Einstein and Stern ([0, p.2): without 
it the energy will depend on the boundaries. Alternatively, if we assume the validity 
of "Kirchhoff's theorem" at high temperatures, we may deduce the zero temperature 
(relative) Casimir energy as follows. 

Let us write the energy per mode k of the radiation field (constrained between 
several disjoint boundaries) as 

E k (T; T c ) = E k (0; T c ) + E' k (T; T c ) (47) 

where E' k is the energy held in the mode without the zpe. For an allowed mode k 

E' k (T;T c )= ^ fc (48) 

exp(pnuj k ) — 1 

(see (fPH) ). Similarly, we denote the analogous quantities for the system with the 
same boundaries far apart as E%°(T) and Ef£°(T), respectively. The relative Casimir 
energy is then 

r c d {T- T c ) = J2 E k (T; T c ) - £ E%°(T). (49) 

k k 



3 The difference in mode densities 8p(u>) ( |34| ) need not be positive definite. 
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Use of ( I^TP ( i.e., separating the thermal energy from the vacuum's) and "Kirchhoff's 
theorem" ( p8|) implies 



hm E^(^^)-E4°°(^) + ^(o) =o. 

T ^°° Ik k J 



(50) 



with £™ l (Q) being the relative Casimir energy at zero temperature (which, of course, 
for the parallel plates geometry, coincides with (|])). Hence 



Now both terms in the square bracket are readily calculable - their temperature 
dependence assures convergence - and thus yield the relative Casimir energy at T = 
without recourse to a (perhaps) objectionable regularization scheme. 

To summarize, the validity, in the classical limit, of "Kirchhoff's theorem" (|28Q 
implies that, in this limit, the relative Casimir entropy is independent of the temper- 
ature T (and obviously, of Planck's constant h), and is thus a geometry dependent 
quantity that characterizes the constraints. Further, its validity allows a novel way 
to evaluate the T = value of the Casimir energy, as we discussed above. 



lim T,E'kiT;T c )-Y,E' k -{T) 

T/TWoo [ k k 



(51) 
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3 Classical Limit for Arbitrary Geometry 



Our discussion of the classical limit of the Casimir effect concentrated so far on 
one concrete example, namely, the parallel plates geometry. However, the intuitive 
physical argument presented in the previous section for the vanishing of the Casimir 
energy in the classical limit is quite general and does not seem to depend in any 
specific manner on the particular geometric details of the parallel plates system. 
Thus, we expect that it should be possible to supplement this argument by an explicit 
calculation of the classical limit of the Casimir effect for a cavity of arbitrary shape. 
In this section we carry out such a calculation, using path integration techniques. 
Some technical details are deferred to Appendices B-D. From this point onward, we 
reset h,c and k B to 1, unless explicitly stated otherwise. 

3.1 Zero temperature analysis 

In order to set up the formalism we start with the zero temperature case. For sim- 
plicity, we concentrate on a real massless scalar field with action 



instead of the electromagnetic field. (Discussion of the electrodynamic Casimir effect 
follows essentially the same lines.) We assume that the field is confined within a static 
cavity of an arbitrary shape and is constrained to vanish on its boundary S, which 
may consist of several disjoint pieces. In addition, we assume that the curvature of 
£ is regular everywhere. 

3.1.1 Sharp boundary conditions over all modes 

Following 0, [7|, H, we parametrize the world- volume swept out by our static £ as 




(52) 




(53) 



where u = (1*1,1*2) are internal coordinates on S. 
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The in-vacuum to out-vacuum transition amplitude Z associated with fl52| ) with 
<f>(x) constrained to vanish on (|53[) is given by 

Z c = i- / V<j> n ^ x(u))] exp^ f d'x (d^) 2 (54) 
with normalization 

Z = jV(j> exp^ | rf 4 x (^0) 2 . (55) 

We chose to normalize ([54]) relative to the unconstrained system, since, after all, we 
are interested only in the difference of the ground state energies. 
Following pi we write 



Here 



IJ S[<j>(t,x{u))] = J Vip(t, u) exp i J dt J^d 2 u^g^(t,u)4>(t,x(u)) . (56) 



d\ 2 / r\ \ 2 / o o\2 
x \ / ox \ /ax ax 



51 \dui) \du 2 j \du-L du 2 j ^ 
is the determinant of the induced metric on E. 

Since E is static, it is convenient to resolve and ip into frequency modes as 

oo oo 

0(t,x)= / -% w ( x ) and V(*,*0 = [ ^ ^ w (x) e~** (58) 

J Zn J In 



-oo 



(with 0* = and similarly for ip). 

Substituting (EBI) in (E3) and using (|5*8"|) we have 



-oo 

(59) 

We may summarize (^) and (|59|) by saying that the ^ act as sources for the 
</>* , which live on E. Integration over these sources constrains all modes (p^ (i.e., low 
frequency modes as well as high frequency modes) to vanish on E. This is, obviously, 
an idealized situation, e.g., if we try to take fl5*2|) more seriously as a simplified version 
of the electrodynamic case. Indeed, it is well known that real metals are transparent 
to radiation of frequency uj > u p , where uj p is the plasma frequency of the metal in 
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question (which typically lies in the UV region). Low frequency modes (i.e., modes 
with uj < u p ), on the other hand, penetrate the surface of the conductor only to within 
a distance of order l/u p (the "skin" depth) and are otherwise excluded from the bulk 
of the conductor. Thus, as is well known, the "perfect conductor" boundary condition 
that all modes w (x) vanish on E, is an effective description (or a consequence) of the 
microscopic interactions of (j)^ and its sources ip w in the medium, where we assume 
uj p — > oo. The plasma frequency u> p serves as an effective UV cutoff frequency for the 
Casimir effect. In the next sub-section we will treat this issue in some detail. One 
of the benefits of doing so, besides gaining better physical understanding, is that the 
formalism we follow will enable us to calculate deviations from a perfect conductor 
behavior in a straightforward manner. 

3.1.2 More realistic interactions at the boundary 

Thus, to make our model more realistic, we should really think of the conducting 
surface £ as having some thickness w of a few 1/uj p . (In doing so we make the tacit 
assumption that local curvature radii everywhere on E are much larger than l/u p .) 
We denote this thickened three dimensional conductor B. Only modes of frequency 
uj < uj p are constrained to decay into B. Modes of frequency uj > uo p will not be 
subjected to such a boundary condition. This physical behavior should be taken into 
account in (|59|). 

Thus, we replace fl52j) by the action 



where, in analogy with electrodynamics, we introduced the dielectric function 

r 1 + 47TXH, xe£ 

e(w,x) = < 

{ 1, x^S 



with susceptibility x(o>), which encodes the polarization effects due to the W — if) t 
interactions in B. 



(60) 
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As discussed in the previous section, the cavity bounded by B has an IR cutoff 
frequency u c , which obviously sets the scale for frequencies that dominate the Casimir 
effect. We assume henceforth that the ratio uj p /u c is large but finite so as to leave a 
wide window open between the two cutoff frequencies. For the parallel plate capacitor 
of Section 1, with plate separation d, we have of course 

u c = ir/d. (62) 

Then the condition uj p /uj c = Upd/n » 1 means that we are allowed to use the perfect 
conductor boundary conditions only for plate separation d in the range 

d » 0.06/im (j^j ■ (63) 



If we decrease the plates' separation d below what is allowed by (p3|), lu c would tend 
to uj p until it would eventually cross it over. Field modes around u c would penetrate 
yet deeper into the boundary, until it would become effectively transparent to these 
modes, rendering the boundary conditions less effective. Thus, we expect that if our 
modified model (as given by fl6"0"|) and (pID) could still be trusted when oj p /uj c < 1, 
then the — l/d 4 behavior of the Casimir force (10) would crossover to a less divergent 



behavior. In fact, assuming for x( u )i we can show that the Casimir force behaves 
like -1/d 3 for d below the range allowed by fl63|) (see also QT2J , p.230). 

Current experiments are performed (albeit not in the parallel plate geometry) 
with d > 0.1/im|19[]. Obviously, it would be interesting to check experimentally 
(when technology would allow performing such experiments) how the Casimir force 
behaves at such much smaller distances (which are yet large enough on the atomic 
scale) . 

In addition to fl63|) , we assume that u c is large enough, so that for u > uj c we may 
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neglect all resonances and dissipative effects encoded in x, and thus writef] 

4ir X (u) = - < 4, (64) 



in direct analogy with the theory of dielectrics. (In reality, this corresponds to fre- 
quencies in the infrared and higher. In the case of the parallel plate capacitor, tuning 
uj c to the infrared, means plate separation d around 1/xm.) 
With (|64]) for x, the action (|60|) is equivalent to 



2 J J 2tt 



uj 2 + V 2 - m 2 (x) 



with 



m 2 (x) 



pi 



(65) 



(66) 



0, x^S. 

4> develops a finite mass u p inside E and is thus expelled from it.0 Moreover, we see 
from ( |65D that in the limit u p — > oo of a perfect conductor, ^(x) tends to vanish on 
E. 

e(uS) — 1 = 47ix(uj) in (|60|) is, of course, the Green's function (at k = 0) of the 
source field ip (already mentioned after Q53p), that lives in B. Thus, ( |5UD is equivalent 
to 



S, 



eff 




- / — / d*x 0: 



0J 



+ I d 6 x 

IB 



+ (Ci + h.c. 



Thus, we replace (j59|) by 



Yl V^Vipu exp iS'e//, 



(67) 



(6* 



w>0 



9 Following elementary dispersion theory, we may take (as the simplest example) x( w ) = 
—jcr/u)(u) + iy). The simple pole at ui = corresponds, of course, to B being a conductor, with con- 
ductivity a and damping coefficient 7 = u)p/4ir<j. Our second assumption is simply that lu c >> 7. 
More generally, we assume that uj c is much bigger than the highest resonance frequency and largest 
damping coefficient appearing in %. Finally, note that the analytic structure of x(w) insures that 
the path integral over (60) is well defined, namely, Ke(iS) < always. 

10 Long ago, Anderson |45| pointed out that the electrodynamical analogs of ( |65| ) and (|6^) may serve 
as a low energy, non-relativistic example of dynamical mass generation of gauge fields, analogous to 
mass generation in the Schwinger model or to the Higgs phenomenon. 
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with 



oo 



W n °* «pi /£/,■'■* ■ (69) 



With (164]) , the quadratic piece in ^ in (|67|) is +| , w | 2 /a;p. In the limit u p — > oo of a 
perfect conductor, +|-0 ti; | 2 /u;? goes away, the thickened conducting surface £> shrinks 
back into E, and ([BSD reduces to (|59|) (as we already commented following (|66|) ). 

Observe that in our discussion of the Casimir effect in Sections 1 and 2, and in this 
section up to fl59|) , we did not make any reference to the fine structure constant a (or 
its analog, associated with the (pip coupling in £>). However, since io 2 v is proportional 
to a, it rightfully sneaks into our discussion following the modified action ( |60|) , or 
equivalently, (|67|) . Given this, we should, in principle, include the (pip interaction (or 
analogously, the photon-electron interactions) also in the bulk part of the action in 



(|67|) , taking into account the full dynamics of the source field ip. This would imply 
using below the renormalized (p (or photon) propagator in the bulk (for such radiative 
corrections in the Casimir effect see, e.g.,|38|]). However, such effects are small and 
we ignore them in this paper. 

Eqs. ( |67D and d68|) capture some aspects of the dynamics of ip (and thus, in a 
sense, of the dynamics of the boundaries where ip lives), rendering our formalism more 
realistic as it relaxes the constraints on higher frequency modes and thus regularizes 
the theory. One aspect of this regularization is that the two dimensional surface 
£ becomes "fuzzy" and swells into a shell B of finite width. In this context we 
should also mention a recent work by Ford and Svaiter which also takes boundary 
dynamics into consideration, but in a different (and presumably complementary) way. 
In the boundary as a whole is quantized - its location serves as a collective heavy 
coordinate which sits in its ground state (at T = 0) and is thus not sharply defined. 
This quantum uncertainty in the boundary location also regularizes the theory since it 
cures an old problem having to do with divergences of the energy density near sharply 
defined boundaries (see, e.g. the work by Candelas and Deutsch (1979) in The 
idea is simple: if the field vanishes on a sharply defined surface, its momentum at all 
points on the surface is unbounded due to the uncertainty principle, rendering the 
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field energy density singular on the surface. If the location of the boundary has zero 
point fluctuations, the field at the ( "fuzzy" ) boundary is not sharply constrained, and 
the field energy density is no longer singular. 

In the spirit of |J [7|, U we would like now to integrate (p out of ( |68|) . To this end 
it is useful to extend the definition of the sources ip u from B into all space as 

r ^(x), x g b 

VUx) - (70) 
[ 0, x^i3. 

In terms of these extended sources, clearly all spatial integrals in ( |67|) involving <p, 
including the (pip interaction term, are carried over all space. Thus, integrating <p 
out is straightforward and produces a factor Z , which cancels the corresponding 
normalization factor in (p3). We end up with 



Zr 



1 



H Vip u exp 



duo 
2^ 



d 6 xd 6 x' ip*(x.) 



x,x'&B 



5® ( 3 




where 



G w (x a ,x 2 ) =< xi| 



x 2 >= 



g-i|a;| |xi-x 2 | 

4:71 |Xl — X 2 | 



(72) 



In order to proceed, we recall at this point that B is essentially the thickened 
two dimensional conducting surface E, with normal width w of a few skin depths 
l/ujp {I ~ 0{1)). Thus, we parametrize B by 

x = x(u, w) , 



where u are the two internal coordinates of E (see (p3|) ), and w is the coordinate in 
the direction normal to E. We are interested, of course, in the limit of large u p , or 
very small width w. Clearly, in this limit ip u (as well as (p w ) do not fluctuate in the 
normal direction. Consequently we may write the -^-quadratic action in flTT| ) as 

2 

(73) 



S<, = 1 j j 'd*u dV C(u) ^ 5(2) ("-"')+(^) ^(u,u') 



24 



where 



MJu,u) 



g(u) G w (x(u),x(u')) Jg(W) 



9W 



e -iH |x(u)-x(u')| 

4/T |x(u) — x(u') 



is the kernel of a symmetric operator 

We see that Z c is expressed entirely in terms of the sources if) u 
readily integrate over the ip^ and obtain 

1 



c>o det 



cxp 



r 
2" 



da; 



trlog l s + Iuj p M u 



(74) 



We may now 



(75) 



Here T (implicit in the product in the first line of ([75])), is a temporal IR cutoff which 
arises due to time translational invariance. The prime to the right of the product 
symbol after the first equality in ( j75|) is to remind us that t/>( w =o) is a real field, and 
thus integration over it yields a square root of the appropriate determinant. This is 
taken care of in the last expression in ([75]), where we used M.-^ = M-u- The Casimir 
energy is thus 



duj 
— 1 
2vr 



trlog l s + lu p M L 



(76) 



Eq. Q7T3P is one of our main results and we will use it to infer the classical (high 
temperature) limit of the Casimir effect in the next sub-section. 

Note that if £ is the union of several disjoint two dimensional surfaces S a (a = 
1,2, • • • N), where S Q has local coordinates u^ a ^ and induced metric g^ a \ then obvi- 
ously ipuj is the sum 
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of contributions from the individual E a 's (in obvious notation). It is clear that in 
such a case, Xi^ is actually anAf xiV operator matrix jCi^f with kernel 

Mff(u<- a \uW) = ^W(uW) < x(u^)| — -|x(u (/3) ) > y/gW(uW). 

(78) 

( |5[] actually discusses the most generic case where the field lives in D + 1 dimensions 
and is constrained to vanish on several manifolds E a of various dimensions, some 
of which may be time dependent.) In particular, for the parallel plate geometry 
discussed in the previous sections (Si is the z = plane and E 2 is the z — d plane), 
a straightforward calculation yields || |7|, || 

with 



K = ^-^-dl-dl + ie. (80) 

If we let uj p tend to infinity in (|76[) we obtain the Casimir energy of the perfect 
conductor E, identical (up to an infinite additive constant proportional to the area 
of E) to the Casimir energy one may calculate directly from ([5^). In particular, 
substituting ([HJ) in ( [75] ) yields (after proper regularization, and up to an obvious 
factor 1/2) the Casimir energy (H). 

However, keeping uj p large but finite, we may calculate corrections to the Casimir 
effect due to imperfect conductor behavior by expanding ( |76[ ) in inverse powers of uj p . 



Such calculations were made in |30[ using Green's function techniques. In Appendix 
C, following |5|, [7j, ^] (and the formalism presented in this section), we briefly sketch 
the calculation of the Casimir force at zero and finite temperatures for the parallel 
plate capacitor using ([75]). We also calculate in Appendix C the leading (in 1/cUp) 
correction to the Casimir effect due to imperfect conductor behavior for the parallel 
plate capacitor, which coincides with the results of |3(| and |^6| . 
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3.2 Finite temperature behavior and the Classical Limit 

The finite temperature case is obtained, in the usual manner, from the zero temper- 
ature case of the previous sub-section, by continuing u to imaginary values, where 
it can take only the discrete Matsubara frequencies u n = 2nn/ (3. Consequently, we 

oo 

replace continuous integrals / du) /2n by discrete sums Y^=-oo (with (3 = 1/T). 

— oo 

Obviously, since u p plays the role of a UV cutoff, we must sum only over |o? n | < uj v . 
We show in Appendix B that the physics of the electron Fermi gas in the metal 
actually warrants u)\ << u p , so that these discrete sums include many terms, and 
also that we may ignore the T dependence of uo p . These facts indicate that our 
cutoff procedure, using (|6l]) and (|64]), is physically consistent also at finite and high 
temperatures. 

Under the analytic continuation uj — > — iu, the oscillatory exponentials in ([72] ) and 
Q73D are replaced by decaying exponentials. Thus for finite temperature we replace 
by 

e -\u n \ |x(u)-x(u')| 



M.M-y/gfr) 4n |x(u) _ x(u , )| («> 

(as well as a similar replacement for M.^ in (fT8|)). We also introduce the obvious 
notational replacement Ai^ — > Xi n - Note that in the Euclidean domain, Xi n is not 
only symmetric, but also real. 

In particular, note that at finite temperature we must replace ([79|) and ( |5UD of the 
parallel plate capacitor by 



1 e ~ dhn \ 1 
1 / 2h r , 



M n = J dhn \ ) — > (82) 



with h n = ^juj\ — &l — d%. 

Finally, we analytically continue (f75| ) to imaginary uj (and also replace T — > j3 in 
(|75|) ) and obtain the partition function as 

z c = n ' 1 



n>0 



det 



1 S + luJ p M n 
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= exp l~ jr trlog [l E + Zw p jM„] 1 . (83) 

L n=— oo ) 

Thus, the Casimir free energy (the analytic continuation of (ff6|)) is 

1 1 °° r - i 

T c = -- log Z c = — Y, tr log 1 S + lUpMn • (84) 
P n=-oo 

3.2.1 The Classical Limit 

Our main interest in this paper is of course the high temperature, classical limit of the 
Casimir effect. Recall from Section 1 (see in particular Fig. 1) that for the parallel 
plate capacitor (and electromagnetic fields), the classical limit is 

Z^ir or 

(namely, fl2"5]), which we repeat here for convenience). As discussed in Section 2, it 
is supposed to become valid for temperatures T >> T c , where T c = hcn/k B d (see 
Ql6|) ) is a temperature scale defined by geometry. Since u p is our UV cutoff, it is 
clear that the relevant temperature range must be far below hu p (which corresponds 
to ~ lO 50 ^), 

k B T c < k B T « huo p , (85) 

which is consistent with the condition oj p /uj c » 1 (Eq. (p5[)) and also in accordance 
with the discussion in Appendix B. 

According to the general discussion in Section 2, we would expect this behavior 
to be generic, and occur for cavities of arbitrary shape. Thus, we expect that there 
should be a temperature scale, which we also denote by T c , which is determined by 
the geometry of the cavity, in analogy with (|TED, and marks the onset of the classical 
limit. 

In the rest of this section we prove this expectaion for cavities whose boundaries 
consist of several disjoint pieces, and also set a bound on T c . Cavities bounded by 
a single connected boundary are in principle more difficult to handle in this respect, 
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for reasons that we shall explain below. In fact, for such cavities, we found logarith- 
mic corrections to (p5f), which seem to be consistent with the results of Balian and 
Duplantier [ |U2| . 

We see from (|81| ) that A^ n (u 1 ,u 2 ) decays exponentially fast, unless 

u n |x( Ul )-x(u 2 )| < 1. (86) 

In particular, it is clear from (|81|) and (the analytic continuation of) (|78|), that for 
n > 1 

Mf(uW,uW)<^WjW— , a ?p, (87) 

where d a p is the minimal distance between the disjoint pieces S Q and E^. 

Thus, if T = uji/2tc is such that u>id min > 1, where d m i n is the minimal d a p, then 
all matrix elements of M.^ (for any pair a ^ (3) will be exponentially small. Below 
we show that these Ai" 13 m& ke exponentially small contributions to the Casimir free 
energy already when 

2uj x d min > 1 . (88) 

The off-diagonal blocks of ( |§2] ) are a clear example of this behavior, for which we 
can write the operatorial inequality 

g — dhn g—Und 



<- Is- (89) 

Indeed, as we show at the end of Appendix C, the high temperature behavior of (|89|) 
offers a very simple explanation of the exponential decay of the Casimir energy e c (t) 
mentioned following (^4j) (which we derive in Appendix A using a laborious mode 
summation), and of its corollary that the bulk number of modes per unit volume is 
independent of d. 

The diagonal blocks of (^2|) obviously decay only as powers of l/u> n . This is in 
fact true in general, as we now show. According to (|86| ) we should consider only 
neighboring points x(u) ~ x(u + a), which means that they must belong to the 
same S Q (assuming finite d m i n , of course). By continuity then |x(u) — x(u + a)| 2 ~ 
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E„„, g „,,<W, where g „„, = . £r is the induced metric on E mentioned above. 
From these considerations and from the representation! 11 ] 

<f(ae) = lim i- e ~W/ e , (90) 

we immediately conclude that 

M n (ui,u 2 ) d( x(ui) - x(u 2 ) ) = ; (91) 

u) n ^oo TTU) n Z7lUJ n Ui — U.2 

Thus, from the elementary formula 



f ,2 .2 yy^i) d Ul -u 2 1 f 2 r— — 

J Z7T |Ul — U2| I J 

we conclude that 



M n . (92) 

Eq. ( |92] ) is an asymptotic formula for u n — > 00, so for l/d min < uj n « uj p we 
might still have to include many subleading corrections to fl92j)) but the point is that 
the diagonal blocks decay only as powers of l/uj n . 

To summarize, when T is high enough so that (|88|) holds, the off-diagonal blocks 
-^n' 3 ( a 7^ /^) become exponentially small (for 11 ^ 0) , and are thus dominated by 
the power-like decaying blocks. Thus, it makes sense to split M. n (n 7^ 0) into its 
diagonal and off-diagonal pieces 

M n = M n , d + M n , od (93) 

(in obvious notations). 

This splitting is not only according to the different asymptotic behaviors of M. nd 
and Xi nt0 d at large u n , but also according to their different physical interpretations: 
Consider, for simplicity, a boundary £ which is the union of just two disjoint con- 
nected pieces £1 and £ 2 . (The following discussion may be generalized immediately 
to an arbitrary set of disjoint boundaries E a .) Assume that £1 and £2 are infinitely 



n Note that ( p0| ) is the first term in the expansion l x '/ e = (X)^Lo e 2 ™^™) 5 (x), which may be 
used to calculate higher order corrections to (|9l|). Thus, the leading correction to ( p2|) is ~ V^/w^. 
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separated at first, and thus do not interact. Clearly, in this case, the only non van- 
ishing blocks of Xi are the diagonal ones Xi 11 and Xi 22 , that encode the interaction 
of the field cf) and the corresponding boundary Si and S 2 , respectively. By definition, 
Xi 11 depends only on the details of Si and Xi 22 depends only on the details of S 2 . In 
addition, note that the matrix elements defining Xi 11 and jCi 22 (see([73D) are mani- 
festly invariant under three dimensional rotations and translations. Thus, the spectra 
of Xi 11 and Xi 22 are invariant under rigid rotations and translations of Si and S 2 , 
respectively. Then, as Si and S 2 are brought closer together, they start interacting, 
and the off-diagonal blocks Ai 12 and Ai 21 , which encode this mutual interaction, are 
no longer vanishing. As S x and S 2 are drawn closer, the off-diagonal blocks clearly 
become larger, while the diagonal blocks remain unchanged. (We could have invoked, 
of course, this argument to split Xi into diagonal and off-diagonal blocks also at 
T = 0.) The parallel plate capacitor, with Ai given by (^), provides a very simple 
example for this discussion. 

Given our system with boundaries Si and S 2 , we may study its behavior under in- 
finitesimal perturbations. We can perturb the system in two independent ways. One 
possibility is to slightly deform any one of the two boundaries (or both), but other- 
wise keep the state of the plates relative to each other unchanged. (An infinitesimal 
self-similar rescaling is an example of such a deformation.) The second possibility 
is to move the boundaries rigidly relative to each other, that is, to make a relative 
translation or rotation, without changing their shapes. According to the discussion 
above, perturbations belonging to the first possibility will induce changes in the diag- 
onal blocks and also in the off-diagonal blocks. Perturbations belonging to the second 
possibility affect only the off-diagonal blocks. The gradients of the free energy (Q) 
(or of the Casimir energy (|76|) at T = 0) with respect to these deformations yield the 
corresponding Casimir forces. Again, this is all very clear for the example of the two 
plate capacitor. 

Computation of the energy associated exclusively with the diagonal blocks of Xi, 
that is, the energy of the constrained system with Si and S 2 infinitely far apart, poses 
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some conceptual difficulties. 12 Consider Si for example. Its interaction with the field 
gives rise to the divergent (cutoff dependent) quantity trlog.M 11 , which we may 
interpret as the energy needed to assemble Ej. This energy includes the divergent 
self-energy of the "charges" that make Si (which is, by definition, a short distance 
effect), and interactions among patches of charges on much larger distance scales. The 
latter contain the Casimir energy we sought. Then a conceptual question arises how 
to define a finite physically sensible Casimir energy, separate it from the self-energy 
and thus extract it from the divergent trlog Xi 11 of Si. As an example, consider Si 
to be a spherical shell of radius R (an example that has been and is being studied 
extensively in the literature |3|, f|, |32], |37|, P}]). One may attempt to 



define its Casimir energy by comparing the shell with radius R with another shell 
with a much larger radius. But as we inflate the shell of radius R to compare it 
with the other shell, each point on the shell recedes from all the other points on the 
shell. In particular, infinitesimally separated points on the sphere will recede from one 
another too. Thus, inflating the sphere affects interactions over all distance scales, 
and in particular, on the scale of the UV cutoff - that is, the self-energy. On top of it 
all, and related to the difficulty of separating the Casimir energy from the self-energy, 
there is the more practical question of how to separate the Casimir force acting on 
our isolated Si from the total surface tension (and other similar effects), in order to 
measure it. This is to be contrasted with what happens when we rigidly move Si and 
S 2 relative to each other. Then all points in Si move relative to all points on S 2 , but 
there is no relative movement of points within the same S Q , i.e., there is no change of 
the self-energy of each of the boundaries. Such moves thus induce clear-cut changes 
in the Casimir energy of the system ||. 

Having split M. into diagonal and off-diagonal blocks according to ([53|) we expand 
the free energy in powers of the off-diagonal part 

trlog l^ + luj p M n = trlog lx + luj p M n ,d 



12 



We adapt here a discussion given in 
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-tr 



2 [lx + lUpMn,d 



od 



lujp 

1 S + lu p Mn,d 



M 



n,od 



+ 



(94) 



Note that the linear term in Ai no d is absent: 



tr 



od. 



0. 



1 S + lu p M n ,d 

This expansion (|94]) is well behaved due to the exponential bound (^) on Ai n;0 d, as 
we demonstrate in Appendix D. 

Since the leading correction due to M. nfi d is quadratic, we see from ( |37D (and from 
Appendix D) that off-diagonal contributions become exponentially small at temper- 
atures T such that 2u; 1 <i m j n > f . This is (B3) alluded to above. Thus, we expect from 



?8) that T c be bounded from below 



Aixdr 



(95) 



For the parallel plate capacitor, for example, with T c given by ([16]), we have l/Aird = 
T c /4vr 2 . 

Given a configuration of disjoint boundaries S Q , an interesting quantity is the cost 
in free energy (or cost in energy, at T = 0) needed to assemble that configuration 
from an initial configuration of the same boundaries, but infinitely far apart. In other 
words, we would like to study the relative Casimir free energy ([45]), which according 
to the discussion above is given by 



= w { trl °S [is + lUpMn] - tr log [l S + lUpMnJ } 



(96) 



and is a well defined quantity, not plagued by infinite self-energy terms that cancel 
between the two terms on the right hand side of (|9~6"D . 

For the parallel plate capacitor, T r c el (which we calculate in Appendix C), is, 
as expected, 1/2 the Casimir free energy (|22"D which was derived in Section 1 for 
electromagnetic fluctuations. Thus, more generally, we may interpret (^) as the 
Casimir free energy (or the Casimir energy, at T = 0) due to field modes that depend 
explicitly on the mutual arrangements of the disjoint boundaries. (Such modes are 
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to be contrasted, for example, with modes that live inside any of the boundaries S a 
which is a closed boundary enclosing a cavity. When this E a is a perfect conductor, 
the latter modes are disconnected from the outer world, and will thus contribute 
exclusively to Ai aa .) 

It is T r c el in ( p6|) that in the classical limit T » T c becomes proportional to T, 
with a proportionality coefficient (minus the relative Casimir entropy <S£ fl44T ) in the 
classical limit) which is determined by geometry, in a similar manner to (^). Indeed, 
from flSTD, (|95|) and from the explicit bounds ( D.l ) and ( p.2|) in Appendix D, we see 
immediately that for T > 1/ (47r<i m j n ) all terms in (|9~6l) with n^O are exponentially 
small relative to the contribution of the zero mode. Thus, in this limit 

T 

— {tr log 



rel 



Is + IcopMo 



trlog l s + lu p M ,d } 



(97) 



with exponentially small corrections. This is the high temperature limit of the relative 
Casimir free energy, which becomes proportional to T. Thus, the relative Casimir 
energy S£ (T) (|27|) (i.e., the Casimir energy due to field modes that depend explicitly 
on the mutual arrangements of the disjoint boundaries) tends to zero exponentially 
fast (at least as fast as e - 2 ^ d min^ anc [ "Kirchhoffs theorem" (||) holds. Thus, 



Sl el = ~\ log [l s + lujAo] - trlog [l s + lu p M 04 ] } 



(9c 



is the asymptotic relative Casimir entropy ([$(]) (i.e., the high temperature limit of the 
Casimir entropy due to field modes that depend explicitly on the mutual arrangements 
of the disjoint boundaries). Expanding fl9S|) formally in inverse powers of lui p we obtain 

1 



S. 



rel 



--tr log 



(at uj p 



+ ■■■ (99) 
oo, of course), 



It will be very interesting to prove the consistency of ( 
and the analog of (0) for scalar field fluctuations. 

In Appendix C we calculate (|9~7D - (|9~9"D for fluctuations of the scalar field in a 
parallel plate geometry. As expected, we reproduce there half the classical limit 
values of the Casimir thermodynamic quantities given in (p5|) for fluctuations of the 
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electromagnetic field, and also write down explicitly the leading correction in l/lu p 
to these quantities. 

We end this subsection with a remark concerning the high temperature behavior 
of the diagonal blocks of A4. (We remind the reader that these diagonal blocks of 
Xi n encode the interactions associated with each surface separately, i.e., when the 
different surfaces are completely decoupled from each other.) For oj n large, M. n is 



dominated by its block diagonal part M. n ,d and approaches the asymptotic form (|92|) 
ls/2a; n (up to power like corrections in l/u n ), and is thus a short distance effect on the 
surface, which must be related to the high temperature limit of the "electromagnetic" 
self-energy of the surface. Thus, trlog Is + lui p M. n tends to 

tr log 



1 + H. 

2u n 



A ^ht)~ Alog it • (100) 



where A is the total area of all surfaces. Thus, JF C , the unsubtracted free energy fl84|) , 
unlike T r c el \ will contain (roughly oj p /uj c ) terms proportional to 

ATlogfeV (101) 

This logarithmic behavior of T c (for the scalar field) is surely akin to the (renormal- 
ized, electromagnetic) high temperature Casimir free energy (|42|). Since uj p is the UV 
cutoff scale, we expect to have renormalization counter terms such as 



Tlog 



to 

LO. n 



that will renormalize ( |101| ) into (fl2p. It will be interesting to see this happen in a 
concrete calculation. 
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4 Concluding Remarks 



This paper gives the classical - the high temperature - limit of the Casimir effect. 
As a start we studied the simple example of the radiation field confined between 
two infinite parallel conducting plates separated by a distance d, and found that the 
Casimir energy vanishes (up to exponentially small corrections) in the classical limit. 

From that observation we abstracted an argument that in the classical limit (de- 
fined to be temperatures such that the Rayleigh- Jeans, i.e., energy equipartition, 
theorem holds) "Kirchhoff's theorem" is valid - i.e., the energy density of the radi- 
ation field is a function of temperature only. This implies that, in this limit, the 
Casimir energy vanishes. We showed that the zero point energy is required to assure 
the validity of "Kirchhoff's theorem" in the classical limit. Alternatively, assuming 
the validity of the theorem allows the evaluation of the Casimir energy at T = 



without recourse to any regularization scheme. We noted, following [12], that these 
results were anticipated by Einstein and Stern in 1913, prior to the formulation of 
modern quantum mechanics and quantum field theory. These authors noted that 
the high temperature expansion of what we now term Bose distribution function is 

(phu « l) 

k B T- -hu + 0(l/T). 



exp((3hu) - 1 2 
Thus we have here a temperature independent term which contributes to the total 
energy. Its cancellation, i.e., the validity of "Kirchhoff's theorem" in our presenta- 
tion, requires a positive zero point energy, +^Tluj. Thus the removal of the zero point 
energy by considering a "normally" ordered Hamiltonian does not eliminate the need 
for zero point energy. The zero point energy is seen, according to Einstein and Stern, 
to be required for the correct classical limit or, in our terminology, for the validity of 
"Kirchhoff's theorem" . This stems here from the form of the Bose distribution rather 
than the Hamiltonian. 

The study of the (Casimir) free energy in the classical limit led to demonstrating 
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that the Casimir entropy, is (in this limit) a temperature independent constant and 
reflects finite volume corrections in statistical physics. Thereby the proportionality of 
the free energy to the temperature T, that is well known [B| is seen to be related to 
the Rayleigh Jeans limit |22|] rather indirectly - it is a consequence of the interesting 
result that the Casimir force in the classical limit is purely entropic - in fact geometric. 
It is thence predicted that a large class of Casimir free energies, in the high T limit, 
will turn out to be proportional to the temperature - reflecting the purely entropic 
origin of the Casimir force in this limit. The crucial point is the existence of one to 
one correspondence between states of the constrained and unconstrained systems. 

We verified the expected classical limit for the Casimir effect due to a massless 
scalar field constrained by a set of disjoint boundaries of arbitrary geometries. This is 
achieved by subtracting off the contributions of the field correlation function that con- 
nect a boundary piece with itself without interaction with any other distinct boundary 
piece. In other words, our normalization is the partition functions with the boundary 
pieces infinitly far apart. This removes contributions that are independent of the 
inter boundary pieces distances. Our analysis allowed us also to include corrections 
to the Casimir effect due to deviations from perfect conductor behavior. 
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Appendix A: Calculational details for Section 1 



In this appendix we evaluate the Casimir energy through a power series expansion 
in t (t = n^r in the notations of Section 1). This scheme requires regularization for 
each term in the expansion. The result is that, within such an expansion, Kirchhoff's 
theorem is exact, i.e., the Casimir energy vanishes to all orders in 1/T in this high 
temperature expansion. This is interpreted as implying that the classical equipartition 
theorem is robust - the classical approximation, once taken, is exact to within power 
series corrections. We hasten to add that this is yet another example of incorrect 
handling of infinities which is further discussed at the end of this appendix. 

Let us return to the expression for the Casimir energy density at finite temper- 
ature, Eqs. (fl"2|) and ([18]), which, in our dimensionless units are summarized by 



(x m = Vx 2 + m 2 ), 

oo 

s c {t) = J xdx 





— coth( — ) + x m coth( — - 
2 K 2t J ^ K 2t 



dm coth( 



2t 



(A.l) 



We now adjunct to each of the above integrals the "standard" [17, TH| cutoff function, 

f(x m ) = exp [-ax m ] , 

assuring thereby the convergence of the integrals. (We are interested in the a — > 
limit.) 



Assuming the validity of this regularization scheme, we may expand ([43], p. 42), 
1 z z 3 2 2n Bo 



cothz = - + - — — H h , , 

z 3 45 (2n)! 



2n ^2n-l 



(|*| <7T) 



(A.2) 



{Bin are Bernoulli numbers). The evaluation of e c (t) reduces to the evaluation of 



terms of the form [10 



/•oo rgP °° /-oo /-oo roo 

g p = / xdx — f(x) + V" / xdx x p f(x) — / dm / xdx x p f(x). (A. 3) 

JO 2 m=l Jm JO Jm 

The cutoff parameter, a, is set equal to zero at the end of the calculations; the result 
is 



d p+1 

g p = lim(-l)f +1 — — 



1 1 
— + 



a 



1 



2a a 2 exp(a) — 1 a 2 
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Noting that (H3j, p. 1105), 



1/ ^ B n y n 



- = y 

exp(y)-l ^ n! 
and that -£>2n+i = (n > 1), we get g p = for p even, and for p odd, 

_ B p+3 
9P ~ (p + 3)(p + 2)" 

Returning to ( |A.1| ), we see that upon substituting ( |A.2| ) only terms with p even occur, 
i.e., s c (T) = 0. This should be compared with the exponential decay of the Casimir 
energy with temperature that is implied by flTSp, and is depicted in Fig. 1. This 
"robustness" of the classical (erroneous) solution, i.e., the nonexistence of corrections 
to this nil result, is worthy of note. It implies that the correct, low T expression, 
(|TE|), obtainable from the correct quantum statistical mechanical approach, can't be 
obtained by analytic means (i.e., power series) from the high T expansion side. 
The case of p = in the above expansion will be recognized as the sum over the 
constrained number of modes, per unit volume, with the unconstrained number de- 
ducted therefrom. That the result is nil constitutes a proof, based on a particular 
regularization scheme, of the assertion (Section 2) that the number of modes, per unit 
volume, is unchanged upon varying the plates' separation, d. 



39 



Appendix B : Consistency of the plasma dielectric function as a cutoff 



We used in Section 3 the "plasmon mode" as a physical regulator (at zero and 
at finite temperature) in our calculations, through the dielectric function (|H]) and 
(|64]) associated with it. Thus, in effect, we are thinking of the "charge sources" ip^ 
as "electrons" in a metallic conductor B. Since the plasmon mode is a real physical 
excitation in dielectrics, we expect it to be a consistent regulator. It is quite instructive 
to show that this is indeed the case, especially in the finite temperature case. (For the 
sake of clarity, we display in this appendix h's and /cb's explicitly in the appropriate 
places.) 

As is well known, 

Anne 2 ,„ . 

«i = — — > (B-i) 



m e 



where n is the electron density. In most metals, where n/10 23 cm~ 3 = 0(1), %uo p is of 
the order of lOeV, corresponding to temperatures of the order of 10 5 °K, well above 
the melting point of any metal. It is also of the same order of magnitude of a typical 
Fermi energy 

*=l> j »r (b.2) 

of the electron gas in the metal. Indeed, from ( [B.l| ) and (|B.2| ) we have 

" 3.7 L^—) , (B.3) 



flj^p \ ^Avogadro / 

which is 0(1) for all metals. 

The electron Fermi gas is practically degenerate even well above the melting points 
of all metals, and thus 

k B T « e F (B.4) 

in all our considerations. Only electrons with energies close to the Fermi energy of 

the gas contribute to transport. The electron Fermi gas and the gas of photons in 

13 Recall that the behavior of the gas is governed by £ = X 3 n - (120, 000°K/T) 3/2 (n/10 23 cm- 3 ), 
where A = J 2-kTi 2 /m e kBT is the thermal wave length. For typical metals n/l0 23 cm~ 3 — 0(1). 



Thus, for T << 120, 000°^, we have £ >> 1, which means that the gas is in a high density and low 
temperature state. 
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the cavity are in thermal equilibrium. Thus, in the present context, and due to ( |B.3| ). 

%U\ = 2TxksT « hujp (B.5) 



as well. The frequency cutoff normally allows many Matsubara modes. It is the 
physics of the electron Fermi gas in the metal that implies that k^T be much smaller 
than huj p . Happily enough we need not worry about the Fermi gas becoming appre- 
ciably non-degenerate. Thus, in particular, we may ignore any T dependence of n, 
and therefore of u p in our analysis. The metal will melt at temperatures much lower 
than that. The classical limit is reached way below such temperatures. 
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Appendix C : Path integral calculation of the Casimir effect due to a 
massless scalar field in the parallel plates geometry 



In this appendix we evaluate the Casimir energy due to a massless scalar field in 
the parallel plate geometry using the path integral formalism of Section 3 (following 
in part [|, [7|, ||]) and analyze the classical limit. We also discuss briefly deviations 
from perfect conductor behavior. 

As in the text, we take the plates to lie at the planes z — 0, d. Thus, the surface 
coordinates are = (x,y,a), (a = 1,2). According to (j79|) and (|76|) the Casimir 
energy is given by0 



1 f du) 

Sc = o / T~ tr log 
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lujp 
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-dh u 



(C.l) 



where according to fl80|) h u 



-uj" 



d 2 -d 2 + ie. 



As in the finite temperature case (Section 3.2), we can split 
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(C.2) 



2K \ u 1 J 2h u V € J 
Tracing over the and inserting a complete set of momentum states (with 
components parallel to the plates) we obtain formally 



Sr. 



A 

2i 



dud 2 p 



(2tt) 3 

where A is the plate's area, and 



log 



h 



2h 



-uj 2 + p 2 . 



lujp 
2h 



-2dh 



(C.3) 



(C.4) 



We now Wick-rotate to imaginary frequency uj — > —iuj. Thus h — > \/uj 2 + p 2 . (This 
is equivalent to the zero temperature limit of the free energy (0)). Thus, we find 

S r 1 7h 2 dh 



A 



2 J 2tt 2 
o 



log 



2h \2h 



-2dh 



(C.5) 



14 Note the slight change of notations compared to Section 1, where £ c denoted the Casimir energy 
density. 
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This expression is obviously divergent. However, we can remove the divergent part 
by subtracting from it (in the spirit of Section 3.2.1) the self-energy of the plates (the 
contribution of the diagonal part M.^ a) , 



cself 

A 



1 



2aJ ^ trl ° g L ls + H-A^ 



h 2 dh ( lu p \ .„, , 
— log 1 + ^ (C.6) 



which is independent of d. The resulting subtracted expression, the relative Casimir 
energy (|7]) 



S r - ST lf S. 
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(C.7) 



is manifestly convergent. This subtraction is of course the most natural regulariza- 
tion within the present formalism, and seems perhaps more transparent than the 
exponential regulator (0) in Section 1. We can integrate ( |C.7|) by parts and obtain 
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1 + 



- 1 



1 + - M - 



(C.8) 



The perfect conductor limit corresponds to u p d oo. Slightly imperfect conductors 
can be treated by taking u p d in ( |C.8|) large but finite. By expanding (|C.8|) 



m inverse 



powers of lu p d, we can calculate corrections to the Casimir energy due to deviations 
from perfect conductor behavior. We obtainF] 



S, 
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lu p d 



lujpd e u — 1 
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ItUpd 



(C.9) 



The leading term in the last equation, — 14 4 0d3 , i.e., the value of the Casimir energy in 
the perfect conductor limit, is precisely one half of the T = electromagnetic Casimir 



15 This expansion is well behaved, since in effect, the integrand is expanded in powers of 
(u/u>pd)e u /(e u — 1) and of (u / u p d) 2 e u / (e u — 1), which do not introduce poles at u = 0, and due to 



the extra factor of l/(e u — 1) in (|C.9D are also well behaved for large u 
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energy (|9|). Also displayed in ( |C9| ) is the leading correction to the Casimir energy 
due to imperfect conductor behavior. 

The resulting (zero temperature) Casimir force, per unit area, is thus 



fc(d,0) 



i ds: el 1 7 h 3 dh 



c 



A dd 2tt 2 J f-t , 2h 



1 f u du 



e 2dh 



^d±J (^ + l) 2 e«-l 
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1-— - + ■•■ . (C.IO) 



480d 4 \ lu p d 

Note that for any value of lu p d, the denominator in the second line in ( |C.10[ ) has a 
simple zero u = (and no additional zeros), and is moreover manifestly positive for 
all u > 0. Thus, the force (C.1U) is attractive. We displayed in the last line of flC-lOp 
the expansion of the Casimir force in inverse powers of luj p d. The leading term in 
that expansion, 

-dsF- (CU1) 

is the Casimir force in the perfect conductor limit, and it is precisely one half of 
the electromagnetic Casimir force (|10"D. The next term, is the leading correction 
due to deviations from perfect conductor behavior |29|0 f30| , [5J. This correction 
is manifestly positive and thus acts to reduce ( C.ll ), which makes sense physically. 
This expression, unlike (|C.11|) , depends on the unknown parameter /, of 0(1). Having 
this parameter in what is in principle a measurable quantity, is the price we have to 
pay for treating the thickened boundary £>, at large u p , approximately as a "skin 
layer" of width l/u p in Section 3 (see the discussion in the first paragraph of Section 
3.1.2 and also the paragraph following ([72])). As a consistency check, we compare the 
leading correction, (n 2 /60d A ) (l/lu p d), to an analysis of the same problem (for the 
electromagnetic field) in Section 3 of |3(| . The latter calculation is based on analysis of 
electromagnetic Green's functions, and does not make a "skin layer" approximation. 
16 There is a slight numerical error in the leading correction given by J29l. 
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We see that this leading correction is of the same form of the leading correction in 
Eq. (3.22) of |30| , and equals one half of that correction provided I = 3/2, which is 
indeed of 0(1). 

Obviously, we cannot use ( |C.10|) to study grossly imperfect conductors, nor can we 
use it to study the behavior of the Casimir force at small plate separations d that fall 
outside the range fl63|), since we derived ( |C.10|) within the "skin layer" approximation. 
In order to address these issues we will have to take the width of the plates and 
the three dimensional plasma of "electrons" that live in them (as well as perhaps 
additional dispersive effects f47|) into account, but we will not carry this analysis 
here. 

We now turn to finite temperature analysis and the classical limit. According to 
(|84l) and (|8^) the Casimir free energy is given by 
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with h r 
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&i — dl. We concentrate now on T ni the contribution of ui n to 



(|C.12|) . In a similar manner to our analysis at T = (see (|C.5|) ), we find 
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(C.14) 



These expressions are manifestly divergent. However, similarly to our T = 
analysis, we can remove the divergent parts by subtracting from (|C.13|) and ( |C.14|) 
(in the spirit of Section 3.2.1) the self free-energy of the plates (the contribution of 
the diagonal part -M ni d), 

jzself x 
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which is independent of d. 



The resulting subtacted expression, the relative Casimir free energy (f45|) , 
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is manifestly convergent, and actually coincides with half its electromagnetic coun- 
terpart ( p2|) as — > oo. 

We now consider the classical limit. At temperatures T such that 2duj\ > 1, that 
is, for T > T c /4tt 2 , we see that (|C.16|) is dominated by its n = term 
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with exponentially small corrections from the other terms. Thus, in the classical limit 



we indeed obtain an asymptotic behavior of T r c el j A consistent with 
T = case, we expand (|C.17| ) in inverse powers of lu p d: 
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IQnd? \ 

The leading term in that expansion, —jq^, is precisely one half of its electromagnetic 
counterpart in Q25D, with half the corresponding electromagnetic Casimir entropy 

C(3M 



as expected. Also displayed in ( |C18| ) is the leading correction to the classical limit 
(relative) Casimir free energy due to imperfect conductor behavior. They act so as 
to diminish the Casimir entropy slightly. 
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Finally, for the Casimir force (per unit area) in the classical limit we find 
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Appendix D : Exponential bounds on the terms in the expansion in Eq. 

(H) 



It is straightforward to prove the bound 
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(D.l) 



on the leading contribution of M. n ,od on the right hand side of (|94|) . For to n » 1 we 
may estimate the right hand side of ( p.l| ) further by using the asymptotic form 
for Xi n ,d- Thus, we find in this limit that 



<u«| LJ " 
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l + luj p M n ,d ' V u p 2uj n/ \/g a {w 
and consequently we may estimate the right hand side of (p.l 



as 



luip 2cc' n , 



(D.2) 



where A a is the area of S a . In a similar manner we can show that the next (i.e., 4th) 
order in M. n ^ in (0) is bounded (when ( P2p applies) by 
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and so on. 
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